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ABSTRACT 

We present the first results of our followup programme of optical imaging and spec- 
troscopy of a deep 12/im survey conducted with the ISO satellite (Clements et al. 
1999). We find that the objects are typically of fairly low redshift (z~0.1-0.2), but 
with a tail that extends to high redshifts. The highest redshift object is a previously 
unknown quasar at z=1.2. The sample of objects for which spectroscopy has been ob- 
tained forms a complete subsample of the original survey with an R band magnitude 
limit of 19.6. We are thus able to use accessible volume methods to determine the 
luminosity function. We find that the luminosity function is close to the low redshift 
12/im luminosity function determined by IRAS (Fang et al., 1998; Rush et al., 1993), 
but with a possible excess at the highest luminosities, which also correspond to the 
highest redshifts. This excess is compatible with the rapid luminosity evolution of 
L oc (1 -|- z)**-^ claimed by Xu (2000) and suggested by Elbaz (1999) on the basis of 
number counts. These conclusions remain tentative because of small number statistics 
at the highest redshifts. 

Almost all of the galaxies in this survey are emission line objects. We use the 
spectra to determine emission line diagnostics for the underlying ionisation source. We 
find that the majority of objects have HII region-like spectra (roughly 2/3) but there 
is a significant fraction (about 1/3) that contain an an AGN. This confirms suspicions 
that deep mid-IR surveys will be a powerful way to find obscured AGN. Part of the 
survey area has been imaged deeply by the INT Wide Field Camera. Nearly all 12/im 
sources in this area, down to the 3(t sensitivity limit of the original ISO survey, have 
good optical identifications, the faintest being at R=:23. The prospects are thus good 
for obtaining complete spectroscopy for this sample. With such data we should be able 
to make a clear measurement of the evolution rate of this population and determine 
the role of obscured AGN in the mid-IR. 

Key words: galaxies:infrared, galaxies:starburst, galaxies: surveys, galaxies:evolution, 
galaxies:active 



1 INTRODUCTION 

The ISO satellite (Kessler et al., 1996) has vastly increased 
our observational capabilities in the mid- and far-IR, and 
has ushered in a new series of deep cosmological surveys at 
these wavelengths (eg. Oliver et al. 2000; Elbaz et al. 1999; 
Kawara et al. 1998; Puget et al., 1999; Clements et al. 1999). 



* Based in part on observations carried out at the European 
Southern Observatory, La Silla, Chile, and in part on observa- 
tions with ISO, an ESA project with instruments funded by ESA 
Member States (especially the PI countries: France, Germany, the 
Netherlands and the United Kingdom) and with the participation 
of ISAS and NASA 



In the mid-IR these surveys reach flux limits more than 
100 times fainter than our previous best surveys with IRAS 
(Hacking et al., 1987), whilst also attaining substantially im- 
proved spatial resolution. They thus represent an important 
new tool for observational cosmology. They are especially 
important since there is increasing evidence that much of 
the star-formation activity in the universe is obscured by 
dust, which re-radiates the energy absorbed at optical and 
UV wavelengths into the mid- and far-IR. While the discov- 
ery of luminous and ultraluminous infrared galaxies (LIGs 
and ULIRGs) in the local universe by IRAS was the first 
hint of the potential importance of obscured emission (see 
eg. Saunders & Mirabel, 1996, and references therein), it 
wasn't until the discovery of the cosmic infrared background 
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(CIB) (Puget et al., 1996) and large numbers of discrete 
sources in SCUBA submm surveys (eg. Bales et al., 1999) 
and ISOPHOT far-IR surveys (eg. Puget et al., 1999) that 
the cosmological significance became clear. It now seems 
that such obscured objects might contribute up to 2/3 of the 
integrated emission in the universe (Gispert et al., 2000) and 
might thus be the sites for the majority of star formation. 
Alternatively, a significant number of the obscured sources 
might be powered by active nuclei, and the sources uncov- 
ered by mid and far-lR surveys might thus tell us much 
about the formation of quasars, and thence the cosmic X- 
ray background (Almaini et al., 1999). 

While much of the emission from obscured sources will 
be emitted in the far-IR, at wavelengths >50/Lim, there are 
substantial advantages in pursuing surveys in the mid-IR, at 
wavelengths between 10 and 20/im. These are primarily to 
do with the limitations of the current generation of instru- 
ments. Since ISO had only a 60cni mirror, the angular reso- 
lution attainable at long, far-IR, wavelengths is quite poor, 
leading to problems with confusion and with identifications 
for foUowup surveys. These difficulties are largely avoided in 
the mid-IR The mid-IR region was also the most sensitive 
on ISO, with improvements up to lOOx over IRAS, while 
the far-IR typically only reached lOx greater depths, partly 
as a result of confusion noise (Puget et al., 1999). Despite 
operating at shorter wavelengths, the mid-IR is still capable 
of detecting obscured star formation regions. For example, 
the Antennae galaxy has optically prominent galactic nuclei, 
but mid-IR images, obtained with ISOCAM (Mirabel et al., 
1998), show that the bulk of the star formation resides be- 
tween these nuclei, in a region of high obscuration. The peak 
of the mid-IR flux also coincides with the peak of both the 
850/u,m emission, as seen by SCUBA, and the peak of the 
molecular emission, as revealed by CO observations (Haas 
et al., 2000, and papers in Hughes & Lowenthal, 2000). 

Mid-IR observations, however, do have some problems. 
These are principally the result of the complex spectral en- 
ergy distributions (SEDs) in the mid-IR band. For starburst 
galaxies, which are likely to be the majority of the objects 
found in a deep mid-IR survey, the SED is dominated by 
the so-called Unidentified Infrared Bands (UIB), which pro- 
duce prominent broad features at wavelengths of 6.2, 7.7, 
8.6, 11.3 and 12.7 /im, and cam vary in strength from ob- 
ject to object. Furthermore, there is also the possibility of 
emission and absorption line contributions, the strongest of 
which will be the Si absorption feature at 9.7/im. The effect 
of these features on cosmological surveys in the mid-IR hcis 
been extensively studied by Xu et al. (1998). They find that 
the effects of these features will be more pronounced the 
narrower the observational passband. While most deep cos- 
mological mid-IR surveys have been conducted in the 15/im 
band, the present work uses the broader 12/im filter. The 
Xu et al. modelling clearly shows that this filter will be less 
susceptible to the effects of the complex SED. 

The current results of deep mid-IR surveys with ISO are 
summarised by Elbaz et al. (1999; but see also Oliver et al., 
2000; Serjeant et al, 1997; Clements et al., 1999; Serjeant 
et al., 2000). Substantial evolution in the number counts, 
beyond no-evolution extrapolations from the local popula- 
tion, has been found (eg. Elbaz et al., 1999). The meaning of 
these counts, though, is still unclear. Elbaz et al. (1998), for 
example, claim that a new population of objects becomes 



significant at flux levels of < ImJy at 15/im, producing a 
distinct bump in the number counts. They suggest that this 
is a new population of objects at z<1.5 contributing a sub- 
stantial amount of the star formation history of the universe. 
Meanwhile, Xu (2000) suggests that a combination of simple 
luminosity evolution at a rate of L oc (1 -|- z^'^ combined 
with the colour effects of the UIB features can produce the 
jump in number counts seen in the data. It is not possible 
to decide between these possibilities using the existing rmm- 
ber count data. The only way to determine the true nature 
and significance of the sources uncovered by deep mid-IR 
surveys is to determine their redshifts, and thus to explic- 
itly measure their luminosity function and evolution. Vari- 
ous models of this population can then be explicitly tested. 
We cannot yet directly obtain redshifts for mid-IR sources 
in the waveband of detection, but must instead use ground- 
based optical and/or near-IR observations. Identifications, 
and spectra, at these complementary wavelengths must thus 
be obtained. 

We here report the first results from identification 
and redshift measurements for a deep 12/im ISO survey 
(Clements et al., 1999) conducted using the ISOCAM instru- 
ment (Cesarsky et al. 1996) on the ISO satellite. The survey 
covers a total area of 0.1 sq. degrees made up of four sep- 
arate subregions separated by several degrees, and reaches 
a 3cr flux limit of ~300/iJy. A total of 148 sources were de- 
tected above the 3(t limit, though a number of these may 
still be fragments of brighter sources. Where optical identi- 
fications in the US Naval Observatory Catalogue (USNO-A) 
are available, we can achieve star-galaxy separation using 
an optical-IR colour-colour method. This relies on the stel- 
lar spectra approximating a single temperature black-body, 
while the galaxies are assumed to have a substantial infrared 
excess. Very few of the galaxies found in the survey were 
previously known. Only one object, which was detected by 
IRAS, had a measured redshift. A full description of this 
12/im survey is given in Clements et al. (1999) (hereinafter 
Paper 1). The goal of the present paper is to present the first 
results of the foUowup programme aimed at this survey. 

The rest of this paper is organised as follows. In the 
next section we describe the sample of objects observed, and 
after that the observations themselves. The results of the 
spectroscopy arc then discussed followed by brief comments 
on the identification of 12/xm sources in two previously blank 
fields. Discussion of the luminosity function of the 12/L(m 
sources and then emission line diagnostics for their power 
sources follows. We then discuss the nature and significance 
of this population before drawing our conclusions. An Ho of 
75kms~^Mpc~'^ and qo = 0.5 are assumed throughout. 



2 THE SAMPLE 

The target list for the spectroscopic foUowup was con- 
structed from the list of > 5(7 non-stellar sources given in 
Paper 1. From these we then selected those objects with an 
optical ID brighter than R=19.6 magnitudes, as cross iden- 
tified in the USNO-A catalogue. This provides a list of 25 
objects. [The object F3_15, was found to be a star in paper 
1 on the basis of its optical/IR colours. This stellar identifi- 
cation was inadvertently missed from Table 2 in that paper. 
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but was included for the analysis work in that paper. It is 
thus not part of the basis sample in the current work.] 

As part of the INT Wide Field Survey (McMahon et 
al., 2000) and the UKIRT Mini-survey (Davies et al., pri- 
vate communication), much of Field 1 of Paper 1 was ob- 
served in the optical/near-IR, providing UBVRIJK images 
and photometry. The identifications and fluxes from USNO- 
A can thus be checked and are found to be broadly cor- 
rect. The deeper optical data also allows us to go fainter in 
our identifications of ISO sources in Field 1. We find that 
all of the > 5cr ISO sources are identified in this field (the 
faintest is Fl_44 which is identified with an R band object 
with R=23.15). We also attempt to obtain identifications for 
fainter 3 - 5a significance ISO sources using the new opti- 
cal data. We find good identifications by eye for all except 
three ISO sources > 3cr within the area covered by the INT 
WFS image. These three unidentified sources are all of low 
signal-to-noiso, (3.5, 3.0 and 3.2 a) and so might be residual 
noise spikes or cosmic hits that have managed to escaped 
our ISOCAM reduction. The reliability of the source detec- 
tion, discussed in paper 1, with only 4% of initially detected 
sources being rejected as spurious cosmic hits, suggests the 
that these may be real sources. We thus conclude that the 
optical imaging, to about R=23.5, which identifies 26 ISO 
sources and fails on these three, has a maximum identifica- 
tion incompleteness of ~ 11%. This bodes well for obtaining 
complete identifications for the 3 - 5a sources in the other 
three fields with deeper optical imaging, which has now been 
obtained (Clements et al., in preparation). 

The basic properties of the sources observed in this 
project are given in Table 1. 



3 OBSERVATIONS 

Most of the observations were carried out at the ESO 3.6m 
telescope at La Silla from Nov. 13 to 15 1999. We used 
the EF0SC2 instrument. This provides both imaging and 
long slit spectroscopic capabilities, using filters and grisms 
(Patat. 1999). The detector is a 2060x2060 thinned Loral 
CCD which has good sensitivity throughout the optical spec- 
trum. The pixel scale is 0.157" /pixel. Most of the spectra 
discussed here were taken using a 2" wide slit, though a few 
used a 1" slit, and Grism 12, which covers the spectral re- 
gion from ~ 6000 to 10000 A at 2.12A/pixel dispersion. A 
few objects were also observed in the blue using Grism 11 
which covers ~3400to 7500A at the same resolution. Wave- 
length calibration was achieved using a He-Ar arc lamp at- 
tached to the spectrograph, while flat fielding used a dome 
fiat. The contribution of sky lines to the spectra was cal- 
culated by fitting and subtracting a cubic spline function 
to the spatial direction for each detector row in the dis- 
persion direction. Standard IRAF tasks were used for these 
procedures. The CCD response function was calibrated out 
of the spectra using observations of the spectroscopic stan- 
dard stars LTT9239 and LTT1020. Spectroscopy integration 
times ranged from 900 to 1800s. Typical observing condi- 
tions were for seeing between 1" and 2" , and with reasonable 
transparency. The conditions, however, were not sufficiently 
stable to be photometric. 

Two objects, F1_0 and F2_0 were observed on a previ- 
ous run, from 24 to 25 December 1997, which was severely 



affected by poor conditions and equipment problems. This 
was also conducted at the ESO 3.6m using EF0SC2. The 
observational setup was broadly similar, though a different 
grism, Grism 1, was used. This covers ~3200 - 10000 A at 
6.3A/pixel. The observations and data reduction for these 
objects followed the same scheme described above. The re- 
duced resolution, however, means that these observations 
cannot be used for emission line diagnostics. 

In addition to the spectroscopic observations, R band 
imaging was also obtained for two objects during the 1999 
run. These were F4_2, and F4_ll, chosen to be bright 12/im 
sources, with strong detections in the ISOCAM observa- 
tions, btit without a bright (R<19.6) optical identification. 
The images were taken as a series of four separate 200s in- 
tegrations each shifted by -30" in both RA and DEC. The 
images were flat flelded using an imaging dome flat, and 
them median stacked with the appropriate shifts taken into 
account. Residuals left after the dome flat were them cor- 
rected be median stacking the images for both objects with 
each other without the shifts corrected, to provide a sec- 
ondary, sky flat image. 



4 RESULTS 

4.1 Spectroscopic Results 

The reduced spectra were all visually inspected and initial 
classifications, as stars, emission line galaxies and, in one 
case, a possible absorption line galaxy, were made on this 
basis. The VEist majority of objects turned out to be emission 
line galaxies. These spectra were further analysed to extract 
redshifts, determine emission line strengths, and to search 
for broad emission lines. Redshifts and derived properties are 
summarised in Table 2. The following method was used to 
extract spectral properties: first the brightest pixel in a given 
wavelength range is assumed to be Ha. Then a composite 
spectrum, including a background flux constant wrt. wave- 
length. Ha, [Nil] and [SII] emission lines, is matched to the 
spectrum. The lines are assumed to have Gaussian shapes, 
with [Nil] and [SII] lines having the same width but different 
from Ha. The strengths of the Ha, [Nil] and [SII] lines arc 
allowed to vary, but the ratios of the [SII] and [Nil] doublets 
are fixed at their physical values (Osterbrock, 1989). The 
assumed redshift is also allowed to vary. The between 
the model and the observed spectrum is minimised using 
a downhiU simplex method (AMOEBA in IDL). Where a 
broad line is suspected, an additional model component con- 
sisting of an extra Ha with larger initial Gaussian width is 
added to see if the is significantly reduced. This is found 
to be the case in only one object. The results of these fits 
are given in table 3, as are the results of the spectroscopy 
more generally. Spectra axe shown in Fig 1 and, for the one 
quasar in the sample, F4_12, in Fig.2. 

Given the redshifts we then calculate the 12^m lumi- 
nosities of the objects in the usual manner. We include a 
K-correction to account for the effects of the complex spec- 
tral energy distribution (SED) in the mid-infrared. The K- 
correction adopted is a simple linear fit to the models of Xu 
et al. (1998). Over the redshift range of interest (0 - 1.2) the 
Xu K-corrections can be approximated as K = 1.7 z, using 
the definition of K-correction from Xu et al. (1998). How- 
ever, this is the mean K-correction and does not take into 
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Object 


RA (2000) 


Dec (2000) 


Fi2 (^Jy) 




Rcdshift (z) 


Log(12;um Luminosity) (Lq) 


Log(Bolometric Luminosity) (L©) 


F1_0 


03 05 36.2 


-09 31 21.6 


12147±200 


17.1 


0.113 


10.1 


11.2 


Fl_3 


03 05 15.1 


-09 31 17.6 


3835±99 


18.6 


0.191 


10.1 


11.2 


Fl_4 


03 05 05.3 


-09 32 07.7 


2150±185 


19.5 


0.249 


10.1 


11.2 


F1.5" 


03 05 24.5 


-09 35 50.7 


1046±102 


20.2 


0.478 


10.6 


11.7 


Fl_7 


03 05 39.5 


-09 31 25.7 


968±241 


19.1 


0.114 


9.0 


10.1 


Fl_9 


03 05 35.9 


-09 31 43.8 


3732±142 


18.6 


0.115 


9.6 


10.7 


Fl.lO 


03 05 06.1 


-09 32 41.5 


1903±185 


17.9 


0.113 


9.3 


10.5 


Fl.ll 


03 05 14.9 


-09 31 03.8 


1324±113 


18.9 


0.290 


10.1 


10.8 


Fl_12 


03 05 28.4 


-09 35 17.5 


842±92 


16.9 


0.461 


10.4 


11.1 


Fl_18 


03 05 36.3 


-09 31 32.2 


2190±171 


18.5^= 


0.113 


9.3 


10.5 


F1.30 


03 05 30.7 


-09 33 13.4 


636±88 


18.8 


0.138 


9.0 


9.7 


F1.34 


03 05 39.0 


-09 38 32.2 


563±115 


18.3 


0.098 


8.6 


9.3 


Fl_48 


03 05 11.0 


-09 33 10.0 


605±119 


19.1 


0.120 


8.8 


9.9 


F2_0 


03 01 06.1 


-10 44 27.5 


10479±260 


13.0 


0.032 


8.8 


9.9 


F2.3 


03 00 37.3 


-10 42 51.0 


2385±100 


18.8 


0.171 


9.8 


10.5 


F2_24 


03 00 35.4 


-10 43 22.3 


814±118 


18.5 


0.119 


9.0 


9.7 


F2_80 


03 00 39.5 


-10 39 43.4 


522±89 


19.5 


0.093 


8.5 


9.6 


F3.3 


03 10 00.9 


-08 39 39.2 


1760±95 


14.6 


M star 






F3.4 


03 10 09.1 


-08 38 57.2 


1682±105 


17.6 


0.107 


9.1 


10.2 


F3.5 


03 09 44.7 


-08 32 55.6 


1374d=108 


17.0 


0.132 


9.3 


10.4 


F3_34 


03 09 46.1 


-08 32 44.0 


639±105 


17.9 


0.576 


10.6 




F4.1 


03 04 03.4 


-10 01 18.0 


3670±149 


17.9 


0.175 


10.0 


11.1 


F4.3 


03 03 55.2 


-09 56 59.9 


3207±84 


16.6 


0.171 


9.9 


10.6 


F4.5 


03 03 47.3 


-09 59 18.2 


2481±215 


19.0 


0.383 


10.7 


11.4 


F4_6 


03 04 06.3 


-09 53 47.6 


3051±197 


16.6 


0.214 


10.1 


11.2 


F4_9 


03 03 56.8 


-09 57 31.3 


1234±95 


17.6 


0.171 


9.5 


10.6 


F4_12 


O:-! ():-! 51.7 


-0!) 58 15. () 


<):i'-!=10(i 


19.1 


1.2 


12.0 


12.7 



Table 1. Properties of 12^m Sources 

Source name is given by catalogue number in appropriate field (Fl to F4), as in Paper 1. Objects already known to be stars from Paper 
1 have been excluded. R band magnitudes in field 1 come from the INT WFC survey data, others come from the USNO-A catalogue. 
Signal-to-noise ratio for the 12/im detection can be estimated from the fiux error. Bolometric luminosity is calculated using the 
conversion from 12/im to bolometric luminosity determined by Spinoglio et al., (1993) and the classifications from emission line 
diagnostics. A 'mixed' class object is treated conservatively as being an AGN. Where emission line diagnostics are unavailable the 
object is treated as a starburst. Notes; a: additional object below nominal R band magnitude limit; b; source outside region covered by 
INT WFC survey, R band magnitude from USNO-A catalogue; c insecure redshift based on absorption feature. All other redshifts 

based on emission lines. 



account the effects of SED differences from one object to 
another. Those differences, though, arc smaller in the ISO 
12/jm band than in any of the other ISO filters used for deep 
ISOCAM surveys (Xu et al. 1998). The uncertainties pro- 
duced by this scatter should bo loss than ~0.2 mag, and thus 
should not have a major effect on the luminosity function 
we derive from these luminosities. 



sample. It is also the only object in the sample characterised 
by a clear broad line spectrum, containing emission linos 
of Mgll and AlII/CIII]. There also appears to be a broad 
emission feature to the blue of Mgll which could be due to 
Fell. The Mgll FWHM is found to be 4300fcms"^ (observed 
frame). The spectrum is shown in Figure 2. 



4.. 1.1 Notes on Specific Objects 

F3_3: This object was classified as a galaxy based on its 
location on the B-R-F12 colour/colour diagram (Paper 1). 
However, spectroscopy clearly shows that the object is an M 
star. The previous misclassification is probably due to the 
large molecular absorption bands in the M-stax spectrum 
suppressing optical emission below what would be expected 
from a simple black-body spectrum. 

Fl_34: The Ha emission line in this object appears to be 
broad, ~1000 kms~^ (cf. 400kms~^ for the other, narrow 
lines) indicating the likely presence of a Seyfert 1 active nu- 
cleus. 

F4_12: At z=1.2 this is the highest redshift object in the 



5 BLANK FIELDS 

Ten sources from Paper 1 arc as yet optically unidentified. At 
the same time as our optical spectroscopy, we thus also ob- 
tained moderately deep images for two 12/im optical blank 
field sources, F4_ll and F4_2. These are shown in Figure 3. 
In both cases we find plausible faint optical identifications 
within 7" of the ISO source position. In the case of F4_ll 
the identification is isolated and unambiguous. In the case 
of F4_2, however, there would appear to be a group of sim- 
ilar faint objects close to the primary identification. There 
is thus the potential that the ISO source may be a com- 
bination of two, possibly interacting, objects, or that the 
identification might be one of the other optical sources in 
the field. 

We can obtain an estimate for their magnitude by com- 
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\J UJCC L 


rv,c<_ibiiii L 


JTLCf Hi VV 


1 IN 111 Ijj VV 


[OllJ VV 


HcK (bro&d) 




Fl_3 


0.191 


104.6 


60.7 


26.4 




HII 


F1.4 


0.249 


41.5 


21.0 


15.0 




HII 


F1.5 


0.478 


50.3 


33.4 


<48 




?, HII* 


F1.7 


0.114 


14.7 


4.9 


1.8 




HII 


F1.9 


0.115 


26.7 


11.1 


10.0 




HII 


Fl.lO 


0.113 


17.0 


6.4 


4.3 




HII 


Fl.ll 


0.290 


21.7 


35.9 


<3.0 




M, HII* 


Fl_12 


0.461 


54.9 


53.0 


30.7 




SII, SII* 


Fl_18 


0.113 


79.1 


50.1 


24.9 




HII 


F1.30 


0.138 


28.1 


20.1 


2.7 




M 


F1.34 


0.098 


9.1 


13.6 


3.4 


26.4 


SI 


Fl_48 


0.120 


61.3 


26.7 


21.2 




HII 


F2_3 


0.171 


13.25 


13.3 


8.8 




SII 


F2_24 


0.119 


17.3 


20.2 


17.2 




SH 


F2_80 


0.093 


56.9 


17.1 


18.1 




HII 


F3.4 


0.107 


49.6 


21.0 


13.2 




HII 


F3.5 


0.132 


21.2 


10.9 


A band 




HII 


f4.1 


0.175 


35.0 


20.9 


10.3 




HII 


F4.3 


0.171 


14.1 


11.9 


5.7 




SII 


F4.5 


0.383 


46.7 


44.8 


9.5 




M, SH* 


F4.6 


0.214 


11.3 


6.6 


1.8 




HII 


F4.9 


0.171 


A band 


A band 




? 




F4_12 


1.2 










SI 



Table 2. Emission Line diagnostics and Classifications 

Emission line equivalent widths and derived classifications for objects where emission line classifications are possible. HII indicates an 
HII region— like spectrum, and thus an object likely to be powered by a starburst. SI and SII indicate a Seyfert 1 (ie. broad line) and 
Seyfert 2 spectrum respectively, and thus the presence of, and probable powering by, an AGN. M indicates a spectrum that may be a 

mixture of HII and Sy2, or may be a LINER - the current limited set of emission lines cannot discriminate between these possibilities. 

'A band' indicates that the relevant lines are ailected by atmospheric absorption by the A band, leading to uncertain equivalent widths. 
In most of these cases classification is not possible. ' ?' indicates that the current measurements or limits to equivalent widths do not 
permit classification. In some cases there arc second classifications indicated by * . These arc classifications which include the H/3 and 
[OIII] fiuxes when they are available. In the case of M and ? classifications these later, more complete, classifications are adopted. 



Figure 1. Observed Spectra 

Spectra for the galaxies observed in this programme. The locations of important features, typically Ha and SII, are indicated, and the 
wavelength range shown is chosen to highlight these features. Spectrophotometric callibration is not fully calibrated so flux values are 
indicative only. The objects F1_0 and F2_0 were observed in poor conditions on an earlier observing run. These two spectra are not flux 
calibrated, merely divided through by a standard star, so the apparently blue continuum slope is not real. 



parison with magnitudes of known objects imaged at the 
same time. We find that F4_ll has R~20 and F4_2 R~22. 

These results arc consistent with the results from the INT 
WF survey data discussed above. 

With identifications obtained for these two sources, 
there remain 8 sources with ISO detections > 5ct that still 
lack optical identifications. Since these are the objects with 
the highest 12/L(m to optical ratios, they are clearly interest- 
ing targets. Their optical faintness suggests that they are at 



high redshifts, and they are thus likely to be highly lumi- 
nous. Further observations will be proposed to obtain iden- 
tifications and redshifts for these, and fainter ISO sources. 



6 THE LUMINOSITY FUNCTION 

The data presented in this paper provides redshift informa- 
tion for a complete subsample of the 12/im survey. We are 
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Figure 2. F4_12 Spectrum 
Spectrum for the 12/im selected quasar discovered in this survey. 




Figure 3. R band images of 12/im optical blank fields. 

F4_2 (left) and F4_ll (right). The images are 31.4" to a side. The nominal position of the ISO source is at the centre of each image. 
The astrometric accuracy of the ISO observations is 6" (2it) (Paper 1), and this is shown by the large circle. A smaller circle highlights 
the likely identification. The F4_2 identification has R~22, F4_ll has R'~20. Sec text for details. Images are of similar depths, reaching 

R~23. The scaling is linear, chosen to highlight relevant features. 



complete to the 5(t flux limit of the survey, and to R=19.6 
magnitudes for the optical identifications and redshift in- 
formation. We can thus use the accessible volume method 
(Avni & Bahcall, 1980) to determine the 12/im luminosity 
function in this sample, whose 12/Lim flux limit is a factor 
of > 20 deeper than existing determinations of the mid-IR 
luminosity function (Hacking et al., 1987). We use the stan- 



dard 1 IN max method to determine the luminosity function, 
with the volume determined from the maximum distance 
at which a source would be detected ie. the accessible vol- 
ume. The luminosity function contribution of each object is 
1 IN max where Nmax is given by the redshift at which an ob- 
ject would drop out of cither the survey data or the followup 
programme, whichever is smaller. 
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Thus: 

$(L)dL = ^l (1) 

i 

where is the number of objects Mpc~^ in the lu- 

minosity range L to L + dL and 

Vi = min{Vopt,i, Vi2,0 (2) 

where Vopt,i is the volume in which the iih object could 
have been detected down to the R=19.6 magnitude limit of 
the optical foUowup, and Vi2,i is the volume in which the 
ith object could have been detected down to the 500/iJy 
flux limit in the original survey. The volumes are calculated 
taking into account the variable flux depths of the original 
12/xm survey (see Paper 1). 

The luminosity function produced from these consider- 
ations is shown in figure 4. As can be seen the luminosity 
function is a fairly close match to existing 12/u,m luminosity 
function determinations based on low redshift IRAS data 
(see eg. Fang et al 1998). The only substantial excursion 
from the local luminosity function occurs at the highest lu- 
minosities. Since these high luminosity bins are also dom- 
inated by the highest redshift objects this might be taken 
as evidence for evolution. Combining the two highest lumi- 
nosity bins to provide better statistics (3 objects) at the 
highest luminosities suggests a, ^ 3a significance difference 
between the current luminosity function and that derived 
from IRAS (Fang et al. 1998; Rush et al., 1993). The mean 
redshift for objects in this bin is 0.7, implying density evo- 
lution at a rate of (1 -I- z)® or luminosity evolution at a 
rate of (1 + z)'' with large uncertainties on these exponents. 
These values are broadly consistent with the evolution in the 
ultraluminous IRAS galaxies suggested by Kim & Sanders 
(1998) and with that suggested by Xu (2000) for the 15^m 
population. However, we as yet have only a small rmmber 
of objects at the highest luminosities. Indeed, the highest 
redshift bin contains only one object, F4_12. This is one of 
only two broad-line objects in the survey as well as being at 
the highest redshift, so might be an unusual object biasing 
the statistics. There may also be additional biases entering 
the data, eg. Malmquist or Eddington bias, since the sam- 
ple size is so small. These issues will be addressed when the 
identification programme is complete and we have a larger 
number of sources. 

These hints at evolution in the high luminosity /redshift 
end of the population can only be confirmed by determin- 
ing redshifts for the rest of the sample. Since the remaining 
galax;ies are optically faint (R>19.6) there is a strong pos- 
sibility that substantially more high redshift galaxies will 
be found, allowing a much better examination of the high 
redshift and high luminosity end of the luminosity function. 

7 EMISSION LINE DIAGNOSTICS 

The nature of the ionising spectrum in a galaxy can be de- 
termined by examining emission line ratios (Baldwin et al., 
1981). Standard ratios include Ha/[NII] and Ha/[SII], which 
can be determined for most of the objects in the current 
sample, as well as [OI]/Ha, [OIII]/H/3 and [OII]/H/3. These 
latter ratios are not available for most of the objects in the 
sample since our spectral range did not extend far enough 



into the blue, or go deep enough to detect the weak [OI]6300 
line. Consideration of the standard emission line diagnostic 
plots shows that we can attempt a classification using just 
the Hq, [Nil] and [SII] lines we have for the majority of 
the objects. Our classification diagram is shown in Figure 5. 
On this basis we find that the majority of our 12/jm sources 
(11/16 where appropriate data is available) have HII region- 
like spectra ie. are powered by starbursts, while 4/16 have 
Seyfert2-like spectra ie. contain and may be powered by ac- 
tive nuclei. Two further objects have broad line Seyfertl-like 
spectra, producing a total of 6/18 objects that contain AGN. 
Three objects are not clearly classified by this technique and 
could be composite objects or LINERs. 

We can extend this classification scheme for the four 
objects at high enough redshift and with sufficient signal- 
to-noise in their spectra to detect H/3 and [OIII]. Of these 
objects only one is fully classified using just Ha, [Nil] and 
[SII]. This object, Fl_12, is classified as a Sy2 by both meth- 
ods. Of the other three objects where we can use H/3 and 
[OIII], Fl_5 was previously unclassified due to a poor [SII] 
detection, and the others, Fl_ll and F4_5, have mixed clas- 
sifications. Once H/3 and [OIII] ratios are included we find 
that Fl_5 and Fl_ll appear to be starbursts (ie. Hll-region- 
like). F4_5 appears to have a confused classification, with 
[0III]/H/3 vs. [SII]/Ha indicating an HII region spectrum, 
and [OIII]/H/3 vs. [NIIj/Ha indicating a Sy2 spectrum. On 
further examination it appears that the [SII] line strength 
may be underestimated due to a sky absorption feature at 
the extreme red of the optical band (~9300A). We thus ten- 
tatively conclude that it is in fact a Sy2. 

The final tally is thus 14 starbursts, 7 AGN (Syl and 
Sy2) , 1 mixed and one unclassifiable. Of the objects where 
we have a firm classification, 1/3 (7/21) arc AGN-type and 
2/3 (14/21) are starbursts. We find no significant differ- 
ence between 12/im luminosities of the two classes, HII and 
Seyfert, though the most luminous objects are from the 
Seyfert class. 

The Rush et al. (1993) sample of local 12/xm selected ob- 
jects had an AGN fraction of 13%, so the larger fraction we 
have found here, 33%, may represent some evolution in the 
population. Alternative explainations, such as different K- 
corrections between AGN and non-AGN containing objects, 
greater prominance in the optical for AGN objects, are also 
possible. These possibilities will be investigated when the 
identification programme for the complete 12/im survey is 
finished. 



8 DISCUSSION 

8.1 Physical Associations 

If we examine the redshift distribution, we find several 
clumps of objects that appear to group in redshift. These 
include an apparent group around z=0.11 in field 1, and 
around z=0.17 in field 4. We can test the reality or other- 
wise of these associations by examining the positions of the 
objects on the sky. In field 1 we do indeed find that the ob- 
jects F1_0, Fl_7 and Fl_9 appear to have a close physical 
association. They form a small group covering a region ~ 1 
arcminute across, corresponding to ~ lOOkpc at z=0.11, us- 
ing our assumed cosmology (see Fig. 6). These objects may 
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Figure 4. The Luminosity Function 

The 12/im luminosity function from Fang et al. 1998 (open squares and diamond for two different corrections to the local density field), 

Rush et al. 1995 (with luminosity divided by a factor of 1.5 to correct for photometric problems from the use of ADDSCAN data 
(Alexander & Aussel, private communication)) and from the present work (triangles joined by a line). Note the similarities between the 

older, IRAS LF determination and the present, ISO determination, and the possible departure at high luminosities which equate to 
high redshifts. Also shown (X) are the last two bins of the new luminosity function combined to provide a better determination of the 

evolution at the highest redshifts. 
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Figure 5. Emission Line Diagnostics 

Starburst-type objects will lie in the lower left corner, AGNs in the upper right. Objects in the lower right corner receive a mixed 
classification, and may be composite objects or possibly LINERs. The broad line objects, F4_12 and Fl_34, are not plotted on this 

diagram. 
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Figure 6. Group of Objects in Field 1 

Circles show the 6 arcsccond 2(j astrometric errors on the ISO source position. Text indicates the number of the source in the field 1 
list and its redshift. Image is 100 by 66 arcseconds. Note also the disturbed morphology of source Fl_9. North is up and east is to the 

left. Scaling is linear, and chosen to highlight relevent features. 



thus be part of a small group or weak cluster. One other 
object in field 1 is at roughly the same redshift, Fl_10 at 
z=0.113. It is much further away on the sky, 7.5 arcminutes, 
corresponding to 750kpc. It might thus be part of some weak 
larger scale structure, but this is difficult to assess with the 
small numbers of redshifts we have. A simple analysis, using 
the luminosity function derived in section 6, suggests that 
the association of Fl_10 with the other three objects is only 
marginally significant, at ~90% level. 

In field 4 we find that F4_9 and F4_3 are close in both 
redshift (difiierence <0.0005) and on the sky (39 arcsecond), 
with a separation corresponding to 94kpc at z=0.17. They 
are thus likely to be a close pair of galaxies. F4_l, in contrast, 
is significantly further away in both redshift (0.004 in z) and 
angular separation (4.8 arcminutes). It is thus unlikely that 
there is a strong physical association between F4_l and the 
other two objects. 

It is thus possible that 5 of our I2/im sources appear in 
two groups, though better statistics are needed to determine 
if the correlation function of 12/im selected objects is any 
different from optically selected, normal galaxies. 

8.2 Bolometric Luminosities 

Spinoglio et al. (1995) found a correlation between 12/jm lu- 
minosity and bolometric luminosity. The basis for this corre- 
lation was argued to bo that the wavelength is a 'pivot point' 
between the optical emission, which rises more slowly than 
linearly with bolometric luminosity, and the mid- to far-IR, 
which rises more rapidly than linearly with bolometric lumi- 
nosity. They also found that the constant of proportionality 
relating 12/im and bolometric luminosity is different for star- 
burst and AGN sources — thus Lm = 14Li2 for starbursts 
and Lioi = 5Li2 for AGNe, where L12 = i'{12fj,m)Lv(12nm) 
in solar units. If we assume that our optical classifications 



of these objects genuinely identify the source of their lu- 
minosity, then we can calculate the bolometric luminosity 
of the sources. These arc given in Table 2. We find that 
the typical bolometric luminosity for our 12//m sources is 
around 10^°'^ Lq. The population of mid-IR sources so far 
uncovered would thus appear to have relatively moderate 
luminosities, and to not represent any substantial new pop- 
ulation of dust-enshrouded starbursts or AGN. Our highest 
luminosity source, F4_12, however, turns out to have a bolo- 
metric luminosity of W^^'"^ Lq, approaching that of Hyper- 
luminous infrared galaxies. Since this is the highest redshift 
object in the sample, and our luminosity function analysis 
seems to indicate strong evolution at high redshifts or high 
luminosities, then one can speculate that this might be the 
first member of a lurking population of obscured high lumi- 
nosity systems that we have yet to uncover. Completion of 
our identification programme can test this idea. 

8.3 The Nature and Evolution of the 12/um 
MilliJansky Population 

Low luminosity 12/im selected sources appear to be similar 
to the local samples. We find that about 1/3 of them con- 
tain AGNe while the remainder are powered by starbursts. 
These fractions are somewhat higher than those discussed 
by Rush et al. (1993). At the highest luminosities there are 
hints of strong evolution being detected, but these are cur- 
rently dominated by small number statistics. With deeper 
observations we will not only be able to confirm the presence 
of evolution but also test claims based upon the 15/im sur- 
veys that a substantial new population of objects becomes 
significant at the ^ImJy flux level (eg. Elbaz et al., 1998; 
Aussel et al., 1999). It should be noted that selection of ob- 
jects at 12/im is much less affected by K-corrections due to 
the complex mid-IR spectral energy distributions of star- 
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forming galaxies and AGN (Xu et al., 1998) than selection 
at 15//m. It has been suggested (Xu, 2000) that the bump in 
number counts at ~ 0.4mJy is entirely due to this K- 
correction effect and not due to any new population. While 
the current survey is working at somewhat brighter fluxes 
we find no clear evidence for such a population. However, if 
such sources axe optically faint, whether intrinsically or as 
a result of their redshift, they may well have been missed 
by the optical selection criterion. Completion of the optical 
identifications and redshift measurements for all the objects 
in the survey is needed to test these claims and is thus a 
high priority. 

8.4 Future Prospects 

The present paper is only the first step in following up the 
galaxies found in this 12/im ISO survey. Wc have here dis- 
cussed the results of redshift measurements on 26 of the 36 
galaxies with fluxes above the 5a flux limit of our survey, 
forming a complete subsample with R<19.6. While we are 
able to obtain interesting conclusions, they are clearly lim- 
ited by small number statistics, and we are as yet unable to 
study objects with a large Fi2/Fopt ratio (eg. F4_2 which has 
Fi2 of 3.6mJy, but an R mag. of ~22). Completing the red- 
shift observations for the 10 > 5a sources so far unobserved 
is clearly a high priority. In field 1, where wc have deeper 
optical data, nearly all of the 12/im sources are identified 
with objects brighter than R=23.2. If this is representative 
of the rest of the sample, then complete identifications and 
spectroscopy should be fairly easy to obtain. Perhaps more 
interesting is the prospect of extending our redshift observa- 
tions to fainter 12/xm fluxes. If we go to a 4a detection limit, 
rather than 5a we would include ~ 45 more sources. All of 
these objects, again, are detected in our deep optical obser- 
vations of Field 1, so such a study is eminently feasible, given 
deeper optical imaging and a moderate amount of time for 
spectroscopy. We will thus be able to substantially increase 
the number of redshifts we have for faint 12/im sources, and 
will be able to place the tentative conclusions of the present 
paper on a much firmer footing. Such observations will make 
an excellent parallel study to the 15/im surveys also under- 
way (eg. Oliver et al. (2000); Elbaz et al. (1999)). 



9 CONCLUSIONS 

We have obtained spectroscopy for a complete R<19.6 sub- 
sample of our 12/xm survey with the ISO satellite. The red- 
shifts of these objects range from 0.032 to 1.2, with de- 
rived bolometric luminosities ranging from 10'^ '^ to lO^^ '^ 
Lq. There is tentative evidence for strong evolution at the 
high luminosity and/or high redshift end of the luminosity 
function. We find that 1/3 of these sources contain an AGN, 
while the rest seem to be powered by starbursts. At this 
stage our conclusions are limited by small number statis- 
tics, but there axe very good prospects for extending the 
identification programme for this survey, and thereby allow- 
ing stronger, more definitive conclusions to be drawn. 
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